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ABSTRACT 

An investigation was made of the magnitude of the bearing and dis- 
tance errors which result from two passive directional acoustic sensors 
obtaining a crossfix on a target in a medium with a constant horizontal 
velocity gradient. Results are tabulated in Appendix III. For the con- 
ditions considered, the errors were found to be small in relation to the 
target range but became appreciable in absolute value at very large tar- 
get ranges. However, with a knowledge of the horizontal velocity gradient, 
sensor separation, sensor axis orientation, apparent target range and 
apparent target bearing the errors could be predicted. Further, by 


proper sensor orientation, the errors could be minimized. 
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1. iittreduet 16n* 

Straight-line propagation of sound is to be expected only if the 
velocity of propagation is constant throughout the medium. Because of 
the heterogeneous nature of the ocean, this condition is usually not met 
and marked deviations from straight line propogation have been observed 
to occur. 

The principle factors that influence the velocity of sound in sea 
water and therefore result in refraction of the sound ray are temperature, 
salinity, and depth. Variations in salinity are of importance near the 
mouths of large rivers where appreciable amounts of fresh water run into 
the sea, in the vicinity of large ocean currents such as the Gulf Stream, 
and in water close to the surface where rain and evaporation will have a 
maximum effect. However, by far the most important factors influencing 
variations in the speed of sound are differences in temperature and depth 
of the points through which the sound ray passes. 

A considerable amount of work has been done investigating the re- 
fraction of sound rays caused by variations of temperature with depth. 
The aim of these studies has been to increase the detection capability 
of sensors by placing them in such a manner that "shadow zones" are either 
eliminated or caused to be as far from the sensors as possible. Further, 
the studies have investigated ways to make optimum use of sound "chan- 
nels" caused by having a layer of water in which the sound velocity is 
lower than that both above and below it. 

However, once a detection has been made by two directional acoustic 
sensors and the resulting fix plotted, very little has been done to de- 
termine by how much this apparent fix is displaced from the actual target 


location due to horizontal temperature gradients. Although these gradients 


are normally small in relation to the vertical ones, there are certain 
areas near ocean currents, such as the Gulf Stream, where they do become 
appreciable and thus could possibly cause significant bearing errors in 
crossfixes. Further, as the range at which sensors can detect a target 
is increased through technical improvements and making use of "channel- 
ing'', bearing errors could become significant. 

One investigation of the fixing errors which could result from hori- 
zontal temperature gradients considered the case where there was a sharp 
vertical boundary between two bodies of water having different sound 
propagation velocities. 6 | The next logical condition to consider is a 
constant velocity gradient. Therefore an investigation of the effect of 
a constant gradient is presented in the ensuing sections of this thesis. 
2. Theoretical Background. 

Sound energy takes the form of disturbances of the pressure and den- 
sity of some medium. Therefore, the basic relationships between impressed 
forces and the resulting changes in pressure and density are applicable 
to the theory of sound propagation in water. However, the solution of 
the differential equations applicable to pressure transmission phenomena 
requires a knowledge of certain initial and boundary conditions. This 
knowledge is almost always lacking for sound transmission in Sea water. 
Fortunately, it is possible to discuss some of the important characteris- 
tics of sound propagation in sea water without considering the frequency 
and wave structure of the sound. For such discussions the wave can be 
thought of as a beam or ray, and the discussion is concerned only with 
the geometry of the situation. 

The ray concept is based on the assumption that sound energy is 


transmitted along curved paths, called rays, which are straight lines in 


all parts of the medium where the velocity of sound is a constant, and 
which curve according to certain rules in parts where the velocity of 
sound is changing. However, in making use of the ray path approach to 
the solution of sound propagation problems it must be kept in mind that 
certain assumptions are implicit in this technique. If the conditions 
imposed by these assumptions do not apply, the use of the method is not 
valid. These assumptions are: 

(a) The direction of the ray must not change much over a distance 
of one wave length; that is, the radius of curvature of the rays must 
be large in relation to the wavelength. 

(b) The velocity of sound must not change appreciably over the dis- 
tance of one wavelength. [7 | 

The simplest continuous velocity gradient which can be considered 
is one which is a constant. This particular gradient has the interesting 
property of resulting in a sound ray path which is the arc of a circle. 
This fact is easily shown by a consideration of the following: 

Assume that the velocity of sound, V, at a point (x, y), is a linear 
function of the distance of the point from the origin in the x direction; 
i.e. V_= - + gx, where bi and oA are the velocities of sound at some 
point P and the origin respectively, and g is the horizontal velocity 
gradient, a constant. Figure 1 illustrates the situation under consider- 
ation. 

In the diagram the velocity gradient is in the horizontal (x) direc- 
tion. The gradient at any point P with coordinates (x, y) is, by 


definition, 





By integrating this expression and specifying that the x distances shall 
be measured from a vertical line through the points where the velocity 
is zero it is possible to avoid the introduction of a constant of inte- 


gration and obtain the expression 


(2-1) = Vv 
& 
Then, recalling Snell's Law 


y 
(2-2) ae 


where 8 is the angle measured from the boundary of the medium (an iso- 
velocity line in this case) to the ray under consideration and K is the 
Snell's Law constant for this ray, equations (2-1) and (2-2) can be com- 


bined by eliminating V to obtain 


(2-3) — = — 


Differentiating this expression yields 
K ; 
(els) dx = - som Si ngOed.0 


But, from elementary analyticel geometry dx = tan ®@ dy (note that in this 
situation angles are measured clockwise from the y axis), thus (2-4) be- 
comes 


(2-5) dy = - cos 04 6 


Now integrating (2-5) and specifying that vertical distances, jy, 
are to be measured from the horizontal line through the point at which 
8 = 0 yields 


(2-6) yez- = Sin 6 


By squaring (2-3) and (2-6), then adding them, the following result 


10 





Figure 1. General Ray Path Geometry 
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Figure 2. Constant Gradient Ray Path Geometry 


ll 


is obtained: 


(2-7) ce? = S* 


This is easily recognized as the equation of a circle of radius a 
The circle is centered at the origin specified in the foregoing paragraphs 
as being the intersection of a horizontal line through the point where 
8 = 0 and a vertical line through the points where V = 0. 

Since the Snell's Law constant, K, at any point is the ratio of the 
velocity of sound at that point to the cosine of the angle of inclination 
of the ray at that point, it is readily seen that the radius of the circle 


V 


is actually SS » 1.e. the ratio of the velocity of sound at that 


g cos @ 


point to the gradient times the cosine of the angle of inclination of the 
ray at that point. Thus, if the velocity and angle of inclination of the 
ray are known at some point, the coordinates and angle of inclination of 
the ray are completely specified for the total path of the ray. 

By making use of the fact that a constant velocity gradient results 
in rays that are arcs of circles some useful relationships can be ob- 
tained for tracing the path of a sound ray through the water. Consider 
a circle of radius R with the points (x15 y,) and (Xo » Yo) as in Figure 
2. By geometry, 


x, = R- R cos 85 = R(1 = cos 8.) 


it 
©) 

1 
©) 
Q 
8) 
w 
a) 

it 


R(1l - cos 8.) 


Therefore 


(2-8) P= Ko Xe R(1 - cos 8.) ~ R(1 - cos 8)) 


R(cos 84 - cos 85) 


Now, for a linear velocity gradient, g, the velocity, Vos at a point, 


(x, Yo) s may be found from the following relationship: 


2 


aa | + £~AX 


where vy is the velocity at the point, (x) > y,)- 


Thus 


(2-9) Ascws~ se Zemeby 


Equating (2-8) and (2-9) yields 


ey 


(2-10) R(cos 8, - cos 6,) = : 


Solving this expression for R, 


eed a) 


Sc al g(cos 0, - cos 85) 


But, from Snell's Law, 





V 
(2=12) cos 04 fas 1 cos 0. 
Vo 2 


Substituting (2-12) into (2-11) yields 


ve Vo - V 
(2-4 3) Res 2 1 = i i 
Vy A Vy 
87, cos 8. - cos 84) 3 Gee 2 i - 1) 
- ime : 9 
7 Vy, - V - cos 8 
g cos 0 (t= 2) g Z 
2 V 
2 
Similarly, 
(2-14) AY = Yo - y,; = R sin 8, - R sin 9, = R(sin 0,-sin 9) 


en 1 
= 2 {2 sin = (6, = 8,) cos 5 &, = @)| 
But, from (2-8) 


Ox 


cos 9, - cos 9, 


3 


Thus, substituting this into (2-14) 


: ASX 2g el f 1 
y * \scosmenm=ncosnCrne [2 sin 7082 ~ 81) cos 20) + a) 


ho 
=? sc in lig + 6.) sin ~—(68, - 8_) 
ria | 2 vigeas 2 


: 1 ik 
2 sin Se: = 8,) cos = (6, * 81) 


CO 
sin > (8, - 85) cos = (6, + @,) 
ax |- as ee ae OR | 


| fin 1 
[-sin 568) + 9.) sin 7694 - 8.) 
Simplifying this expression by means of a trigonometric identity yields 


Se 


(72155) y = Ax cot ( 5 =) 


Equation (2-15) is the basic equation used in the model described in 
section 3 which is used to trace the sound rays through the water. 

It should be noted that the minimum radius of the circle, as compu- 
ted from equation (2-13), describing the path of a sound ray is obtained 
when cos 9. and the gradient, g, are both a maximum. This means that for 
the model described in section 3 the minimum radius will occur when @ = 0° 


and the velocity gradient is 0.0005 meter/sec per meter. This radius will 


1500 
(.0005) (1) 


signal (approximately the frequency of a low frequency noise source on a 


be = 3x 10° meters. Since the wavelength of a 200 cycle/sec 
medium speed submarine [4] is about 7.5 meters it is readily seen that 
the conditions of assumption (a), that the radius of curvature must be 
large in relation to the wavelength, are met and the use of the ray path 
method is valid. Further, the conditions of assumption (b), that the 
velocity of sound does not change appreciably over one wavelength, are 


fulfilled since the velocity only changes by 7.3 x _.0005 = 2.5x Wer = 


1500 


14 


.00025% in a distance of one wavelength. 
3. The Model. 

In modeling the problem under investigation the following assumptions 
were made: 

a. A detection had been made by two directional acoustic sensors and 
the apparent fix plotted. 

b. There was no change in horizontal gradient with depth, therefore 
the amount of bending of the ray was the same at any depth and the prob- 
lem could be reduced to one of finding the ray path in a horizontal plane 
only. (It is realized that the horizontal gradient does, in actuality, 
change slightly with depth, thus all results are approximations. How- 
ever, to treat the problem in three dimensions would require an amount of 
time exceeding that available for the conduct of the present investiga-~ 
tion.) 

ce. The velocity gradient was a constant in one direction and zero 
in a direction perpendicular to this one. This allowed the use of the 
circle approximation to the path of the sound rays. 

d. The directional acoustic sensor had no bearing error. To add 
this to the problem would obscure the effect of the temperature gradient. 

Figure 3 is an illustration of the situation studied. An apparent 
target was located at XAPP, YAPP, the apparent fix obtained from sensors 
A and B located at Xl, Yl, and XXl, YYl respectively. The bearings of 
the apparent target from sensors A and B were BRGA and BRGB respectively. 
To relate the fix to some common reference point, an "origin" was chosen 
at the midpoint of a line joining the two sensors. The bearing and range 
of the apparent target from this point were BRGP and RAPP. However, due 


to the horizontal velocity gradient, some bending of the sound rays 


fs: 


XAPP, YAPP BRGT 


ERROR 


XACT, YACT 
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Figure 3. Computer Program Geometry 
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Figure 4. Computer Program Geometry (Sensor Axis Rotated) 
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occurred and the actual target was located at XACT, YACT. The bearing 
and range of the actual target from the origin were BRGG and RACT. The 
distance from the apparent target to the actual target was the ERROR and 
the bearing of the actual target from the apparent target was called BRGT. 
RNGERR was the difference between the actual and apparent ranges. 

If the sensors had always been located along an iSovelocity line, 
Figure 3, with the bearings and distances shown, would be sufficient to 
illustrate the geometry of the problem. However, the inclination, by an 
angle OMEGA, of an axis through the two sensors with respect to an iso- 
velocity line resulted in a need to define angles with respect to an 
isovelocity line. By selecting the x and y axes perpendicular to and 
along an isovelocity line respectively it was possible to make the best 
use of the fact that the velocity gradient was a constant in one direc- 
tion and zero in the other. These angles measured with respect to the 
lsovelocity line were assigned Greek names corresponding to the letter 
appearing after the prefix BRG; 1.e. BRGA, the bearing with respect to an 
axis through the two sensors, became ALPHA, the bearing with respect to 
an isovelocity line. The Greek letter bearings equal OMEGA plus the 
appropriate BRG. Figure 4 illustrates the problem when the sensor axis 
is rotated with respect to the isovelocity line. 

Making use of the relationships developed in section 2 relating the 
sound ray path to the arc of a circle, a computer program was written to 
trace the actual path of the sound ray through the water. Thus, given 
initial bearings of the apparent target from two Sensors, it calculated 
the actual target position from the point of intersection of the actual 
ray paths. The program is included in Appendix I. 


The "mechanics" of the program are illustrated in Figure 5, a flow 
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START 
READ INPUT DATA 


COMPUTE 
SENSOR POSITION 


PRINT 
INPUT DATA, OUTPUT TITLES 
SET APPARENT TGT. POSIT. 


SET STARTING VELOCITY 
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COMPUTE INTERSECTION PT. 
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COMPUTE ERRORS 
AND 
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INCREMENT BRGP 
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BRGP = Apparent Target Bearing 

RAPP = Apparent Target Range 

RAPPMAX = Maximum RAPP to be Considered 
SS = Sensor Separation 






BRGP < 350° 


Figure 5. Flow Chart for Program 


i) 


chart for the program. For the reader interested in greater detail a more 
detailed block diagram may be found in Appendix II. The Roman numerals 

in the blocks of Appendix II correspond to the numerals in the left hand 
margin of the program in Appendix I. 

For a given velocity gradient and sensor separation an apparent tar- 
get position was assumed. From this position the apparent bearings and 
distance of the target from the origin and the sensors was calculated. 
Then the actual position of the target was determined by calculating the 
actual ray paths through the water and computing their point of intersec- 
tion. Finally, the errors were calculated. A new apparent target posi- 
tion was then assumed and the calculations repeated for this new position. 
Apparent target positions were assumed for each 10 degrees of bearing 
from an axis through the sensors (except for the O and 180 degree posi- 
tions where no fix is possible) and for various apparent ranges to a 
maximum range of five times the sensor separation. Further, as stated 
earlier, the angle of inclination of the axis through the sensors to the 
isovelocity line was varied in 10 degree increments from 0 to 90 degrees 
for selected sensor separation and apparent range combinations. 

4. Results. 

To evolve a means for estimating the crossfix errors resulting from 
horizontal temperature gradients, the effect of the following parameters 
on the errors was studied using the model described in section 3: 

(a) target bearing 

(b) magnitude of velocity gradient 

(c) sensor separation 

(d) target apparent range 


(e) sensor axis orientation 


20 


The range of values over which the parameters were varied was as 
follows: 

(a) target bearing (BRGP): 10° increments from 10° to 350° relative 
to the sensor axis except for 180° where no fix was possible due to target 
alignment with the sensor axis. 

b) velocity gradient (VGRAD): .0001 to .0005 meter/second per meter. 
These were found to be the approximate maximum velocity gradients which 
exist over long (about 100 km.) distances in the deep sound channel and 
over short (about 40 km.) distances near the surface at the Gulf Stream. 
[5] 

(c) sensor separation (SS): 1 km. to 400 km. These were considered 
to be the minimum and maximum feasible for obtaining a crossfix. 

(d) target apparent range (RAPP): 1 km. to 2,000 km. but never less 
than .1 nor more than 5 times the sensor separation since beyond these 
values the target bearings from the sensors become almost equal. 

(e) sensor axis orientation (OMEGA): sensor axis 0° to 90° from the 
isovelocity line in 10° increments. 

As might be expected, the errors were found to be small in relation to 
the target distance from the sensors but to be appreciable in absolute 
value for targets at long ranges from the sensors or in areas where a 
large velocity gradient existed. It was found that the errors were pre- 
dictable and that they could be minimized by proper sensor orientation. 
The results are summarized in sections 4.1 through 4.8 and depicted 
graphically in Figures 6 through 20. Although all of the computer out- 
put data used in plotting the graphs for the figures is not included, 
samples of the data which illustrate particular effects are included in 


Appendix III. Also, it should be noted that the abbreviations used in 


ZA 


the figures are defined at the beginning of the computer program in Appen- 
dix 1. 
4.1 Effect of Target Bearing. 

Figure 6 illustrates the influence of target bearing (BRGP) on the 
magnitude of the crossfix errors. As the figure and the data in Tables 
1, 2 and 3 show, for a given velocity gradient condition the error was 
very much influenced by target bearing from the sensor axis when the sen-~ 
sor axis was oriented along an isovelocity line. However, when the sen- 
sor axis was at 90° to the isovelocity line, i.e. along the velocity 
gradient line where OMEGA = 90, target bearing had almost no effect on the 
error. The influence of sensor axis orientation is discussed further in 
section 4.5. For the case where the sensor axis was along an isovelocity 
line the errors were seen to increase rapidly as the target bearings ap- 
proached the sensor axis. As discussed in section 2, this was due to the 
fact that the radius of curvature of the sound ray approaches its minimum 
value as the angle which the ray makes with the isovelocity line approaches 
zero. 

4.2 Effect of Velocity Gradient. 

Figure 6 also illustrates the influence of velocity gradient (VGRAD) 
on the magnitude of the cross-fix errors. As the figure and the data in 
Tables 1, 2 and 3 show, the errors were a linear function of the velocity 
gradient, i.e. the velocity gradient of .0005 m/sec per meter resulted in 
errors 5 times larger than those for a gradient of .0001. Although Pig- 
ure 6 depicts the particular case of a 2 km. sensor separation and 10 km. 
target range, the effect which it illustrates was found to be the same 
for all sensor separation and target range combinations examined. The 


plot for the 90° rotation of the sensor axis from an isovelocity line 
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illustrates the linear relationship between velocity gradient and cross- 
fix errors even more clearly than the plot for the 0° rotation of the 
sensor axis. 

4.3 Effect of Sensor Separation. 

The effect of sensor separation (SS) on crossfix errors is not as 
simply described as the effect of velocity and target bearing. As shown 
in Figures 7 and 8 and by the data in Tables 4, 5, 6, 7, and 8 it was 
found that for any given target range, sensor separation had much less 
effect on the errors than target bearing. Also, the minimum error due to 
target bearing generally increased and the maximum error generally de- 
creased as the sensor separation was increased. However, when the sensor 
separation was twice the apparent target range the maximum error due to 
target bearing was at a minimum value and the minimum error was at a max- 
imum value, thus there was the least dependence of error on target bear- 
ing. When the sensor separation was increased to values more than twice 
the apparent target range the minimum error continued to increase and the 
maximum error now increased also. The effect of sensor separation on 
crossfix errors when the sensor separation was greater than twice the 
target range is discussed further in section 4.8. 

An examination of the data in the figures and the tables revealed 
no simple relationship between the magnitude of the error and the sensor 
separation per se. However, as will be discussed in section 4.6, it was 
possible to formulate an estimating relationship for errors based on the 
ratio of target range to sensor separation. Also, the large variation in 
error for the case where the sensors were 10 km. apart and oriented 90° 


to the isovelocity line will be discussed in section 4.8. 


24 





? 
eeetedates So ~ ae 


eauer i 








OMEGA = 0 


SYaLYN = WOU 





180 270 360 BRGP 


9 


Reference: 
Tables 4,5,6 


Effect of Sensor Separation 


7. 


Figure 


ve 




















Ste 


zy 

} 

Je = 
HI = | 
i 





SYaLENOTDI - woud 


OMEGA = 0 








360 BRGP 


270 


0 
OMEGA = 90 


18 


90 


Reference: 
Tables 7,8 


= .0001 


SS = 100 kn. 
we om mem CS m 400 kn. 
RAPP = 500 km., VGRAD 


Effect of Sensor Separation 


Figure 8. 


26 


4.4 Effect of Target Range. 


Figures 9 and 10 show the dependence of crossfix errors on target 
range (RAPP). Samples of the data used to plot the graphs may be found 
in Tables 4, 8-12, and 13. Figure 9a shows the absolute error as a func- 
tion of target bearing for three target ranges and a sensor separation of 
1 km. Figure 9b shows the relative error (the ratio of the error to the 
target range) for the same target and sensor positions. Figure 10 shows 
the relative error as a function of target bearing for different sensor 
separation and velocity gradient conditions. 

In examining the data it may be noted that there was no simple 
linear relationship which related the absolute error at one target range 
to another, thus interpolation is difficult. However, by plotting rela- 
tive error versus target bearing, aS in Figure 9b, it may be seen that 
the relative error was a linear function of target range. 

As Figure 10 illustrates, when the sensors were on an isovelocity 
line, this relationship was found to be generally true for extreme sen- 
sor separations and target ranges also. However, for target ranges less 
than one-half the sensor separation and for target bearings less than 45° 
from the sensor axis the relationship did not hold. As was mentioned in 
section 4.3, the case when the target was less than one-half the sensor 
separation will be discussed in section 4.8. At target bearings less than 
45° from the sensor axis there was less than a 10° difference between the 
target bearings from the two sensors, thus it is unlikely that a "fix" 
would be considered from these two bearings due to the large errors which 
would be induced by small differences in the measurement of either bearing. 
When the sensor axis was at 90° to the isovelocity line, the relative er- 


ror waS a linear function of target range for all target bearings, but as 
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in the case when the sensor axis was along an isovelocity line, when the 


target range was less than one-half the sensor separation the linear re- 
lationship did not exist. 
4.5 Effect of Sensor Axis Orientation. 

As was noted in Figures 7, 8 and 10, the orientation (OMEGA) of the 
sensor axis with respect to the velocity gradient line had a very large 
effect on the magnitude of the errors. In general, errors were minimized 
when the sensor axis was aligned with the velocity gradient line. Fur- 
ther, the errors were practically constant for all target bearings when 
the sensor axis was aligned with the velocity gradient line. The only 
exceptions to this occurred when the target range was less than one~half 
the sensor separation and this case is discussed in section 4.8. 

Figures 11, 12 and 13 illustrate the effect of rotation of the sen- 
sor axis from alignment with an isovelocity line (OMEGA equals 0°) to 
alignment with the velocity gradient line (OMEGA equals 90). The data 
used to construct the graphs may be found in Tables 14 through 23 of Ap- 
pendix III. 

Figure 11 is a histogram showing the maximum and minimum errors, de- 
pending upon target bearing, for targets at an apparent range of 6 km. 
from sensors 2 km. apart. As the graph shows, by rotating the sensor 
axis from 0° to 90° the maximum errors were reduced by a factor of 10 to 
a value approximately equal to the minimum error. Figures 12a and 12b 
further illustrate the changes in the errors brought about by rotating 
the sensor axis. 

Figure 13 depicts the actual and apparent target locations when the 
sensor axis was aligned with an isovelocity line and when the sensor axis 
was aligned with the velocity gradient line. It should be noted that the 


graph is not to an absolute scale, since on such a scale even the maximum 
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error of 68 m. in 6000 m. could not be seen. However, the relative magni- 
tudes of the errors are correct. As Figure 13 and the data in Tables 14 
through 23 show, the actual target was closer than the apparent target 
when the sensors were in a lower velocity medium than the target and fur- 
ther away when the sensors were in a higher velocity medium. When the 
sensors were oriented along the velocity gradient line, however, this re- 
lationship was reversed. 

4.6 Comparison at Constant Reduced Range. 

Reduced range (REDRNG) was defined to be the ratio of the apparent 
target range to the sensor separation. This term was introduced to facil- 
itate estimation of errors for any sensor separation and target range when 
the errors at a different sensor Separation and target range are known. 
Figure 14 illustrates the linear relationship between relative errors at 
one sensor separation and those at another when compared at the Same re- 
duced range. 

Figure 15 depicts relative error as a function of reduced range for 
various target bearings and a velocity gradient of .0005 m/sec per meter. 
As an example of the usefulness of such a plot, consider the following: 

For a target located at a range of 5 km. on a 30° bearing from the 
midpoint of a line connecting two sensors 1 km. apart, Figure 15 shows a 
relative error of about .003 (REDRNG = 5., BRGP = 30). For a target lo- 
cated at a range of 500 km. on a 30° bearing from the midpoint of a line 
connecting two sensors 100 km. apart (at a reduced range of 5.) a predic- 
ted relative error would therefore be .3 (.003 times 100) for a gradient 
of .0005. Thus, using the fact (from section 4.2) that relative errors 
vary linearly with velocity gradient, for a gradient of .0001 a predicted 


relative error would be .06. The absolute error would then be estimated 


34 





OF Oo 

| | 
Os Sound Velocity Increasing ~——_—-» \ 

| Note: RAPP not to scale 
Ay! | « )-XAPP, YAPP 4 

~  @ XACT, YACT (OMEGA = 0) 
x XACT, YACT (OMEGA = 90) 
Gi \ y 


. a 3 y 
‘° | ss a Za 
Qs | 
. S 
SS O 
\ >” 
< = 


aes SS = 2 km., RAPP = 6 kn. 
VGRAD = .0005 
References: 


Tables 14 & 23 


Figure 13. Actual Target Locations 


a5 








qa 


0 


OMEGA = 





360 BRGP 


90 


OMEGA = 


5.0 


REDRNG = 


SS = 20 kn. 


VGRAD = .00025 


km 
km 
km 
km 


—— - —-ss = 10 
5 
2 
1 


— — —SSs = 


55 = 


Effect of Sensor Separation at 


Constant Reduced Range 


Figure 14. 


36 








-4- . —_> ‘ -4 
= —_ * Le ‘2 
} 
os 
: = ees | ee 
i | } 1 
{ 

















' } ! : 
eee ee a ee et te eh ee — 
‘ ! ‘ ! 





REDRNG 


= 1 kn. 
VGRAD = .0005 


Ss 


0006 
nt Oo oO 
AIAN MS 
La’ tad nn A 
oa Go; 
Aon st 
ma NNN 
nw A wn A. 
O@Ggoao 
™ OWN Ee 
mae At AN 
tay ~ A vn wv 
ae 0 eo 
aN OO CO 
nonutw ow 
A, Ar AY AY AY 
Oogodgd ®& 
Sem e eS 
mmaomnmMa 
' i 
lj 
| | 
‘ ’ 
9 
| 


Figure 15. 


Relative Errors at Constant Reduced Range 


oF 


to be 30 km. (.06 times 500). This compares favorably with the actual 
error of 28,151 m. found in Table 7 for this situation. 

Over a wide range of sensor separations, target ranges and velocity 
gradients considered the error estimates obtained in this manner were al- 
ways within 14% of the value computed by the model described in section 3. 
The unusual behavior of the curves in Figure 15 below reduced ranges of .5 
will be discussed in section 4.8. 

4.7 Effect on Bearing Error. 

For targets located such that there was at least a 10° bearing differ- 
ence between the bearings from the two sensors, the maximum error in bear- 
ing from the origin on the sensor axis was 2.2° (Table 13). Except for 
the extreme senSor separations and target ranges these bearing errors were 
less than 1, thus a graphic portrayal of these would not be meaningful. 
However, it was noted that the bearing of the actual target from the appar- 
ent target position was a function of target bearing from the origin and 
was independent of velocity gradient, sensor separation and apparent tar- 
get range except when the range was less than the sensor Separation. This 
case will be considered in section 4.8. 

Figure 16 illustrates the relationship between the bearing of the 
actual target from the apparent target (BRGT) and the bearing of the ap- 
parent target from the origin (BRGP). Four figures have been drawn to 
illustrate the change in this relationship as the SenSor axis was rotated 
from alignment with an isovelocity line to alignment with the velocity 
gradient line. As the figures show, the actual target was always toward 
the lower velocity medium from the apparent target position except for a 


sensor axis orientation along the velocity gradient line. 
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4.8 Effects When Target Range is Less Than Sensor Separation. 


In the foregoing sections numerous references have been made to the 
fact that when the target range was less than the sensor separation, and 
particularly when the target range was less than one-half the sensor separa- 
tion, many of the estimating relationships were no longer valid. 

As may be seen from Figures 7 and 15, when the sensor axis was aligned 
with an isovelocity line, the effect of target bearing on the distance er- 
rors was the same for reduced range values less than 0.5 as it was for 
those above 0.5. However, when the sensor axis was rotated to alignment 
with the velocity gradient line, instead of the errors being almost con- 
stant the errors were greatest at those target bearings where they had been 
the least when the sensor axis was along an isovelocity line. This may be 
seen in Figure 7 and by a comparison of the data in Table 25 with that in 
Table 26. 

Figures 17a and 17b illustrate the change in the absolute errors as 
the sensor axis was rotated. As may be seen in Figure 17a, the absolute 
errors decreased as the axis was rotated but when the reduced range equaled 
0.5, the variation in error due to target bearing remained constant. This 
is contrasted with the behavior, as indicated by Figures 11 and 17b, of 
the error for values of reduced range both above and below 0.5 where the 
variation in the error as well as its absolute value decreased when the 
sensor axis was rotated to alignment with the velocity gradient line. 

Figure 18 depicts absolute distance errors as a function of reduced 
range for reduced ranges from 0.1 to 1. As the figure shows, the maximum 
errors increased in value as the reduced range increased or decreased from 
a value of 0.5 when the sensor axis was along an isovelocity line but de- 


creased uniformly as the reduced range decreased when the sensor axis was 
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along the velocity gradient line. 

Figures 19 and 20 illustrate the effect of apparent target bearing 
from the origin (BRGP) on the bearing of the actual target from the appar- 
ent target (BRGT). The figures also illustrate how the bearing correction 
varied with the sensor axis orientation. A comparison of these figures 
and the data in Tables 24-28 with Figure 16 and the data in Table 8 reveals 
that as the reduced range was lowered from a value greater than one (Fig- 
ure 16) through 0.5 (Figure 17) to 0.3 (Figure 18) the relationship between 
the bearing from the apparent to the actual target and the bearing from the 
origin to the apparent target reversed. 

5. Conclusions and Recommendations. 

As a result of this investigation the following conclusions may be 
drawn regarding the crossfix errors of two directional acoustic sensors 
in a medium with a constant horizontal velocity gradient: 

a. Although the distance errors are small relative to the target 
range, they can become appreciable in absolute value at very large target 
ranges. 

b. For gradients within the values studied, the distance errors are 
a linear function of the velocity gradient. 

ce. The distance errors are sensitive to the ratio of the apparent 
target range to the sensor separation. If two directional sensors are 
used to obtain crossfixes and it is known from past experience that a 
large proportion of their detections occur at a certain range, they can 
be placed so as to minimize the crossf1x error. 

d. The distance errors are sensitive to target bearing and sensor 
axis orientation. By placing the sensors along the velocity gradient line, 


distance errors due to target bearing can be minimized. 


ay 


e. The distance errors are predictable, 1.e. the errors can be esti- 
mated from a knowledge of the horizontal velocity gradient, the sensor 
separation, the sensor axis orientation, the apparent target range and the 
apparent target bearing. 

f. The bearing from the apparent to the actual target is independent 
of velocity gradient, sensor separation, and apparent target range. It is 
only a function of target bearing and sensor axis orientation, thus it is 
also predictable. 

Since the model used to formulate these conclusions did not take into 
account the variations in the path of the sound ray due to vertical velocity 
gradients or horizontal gradients other than linear ones, the investiga- 
tion of errors using a more sophisticated three-dimensional model is recom- 
mended. Such an investigation should, of course, determine the sensitivity 
of the errors to bottom conditions such as slope, since a sloping bottom 
will also result in deaieeaian of the sound ray from its original course. 
Further, since some directional sensor systems depend upon the time phase 
relationship between a pulse received by two sensors, the errors caused by 
the actual paths being of different lengths than the apparent paths may be 
entirely different from those calculated by the model used in this investi- 


gation. Therefore, an investigation of this problem should be considered. 
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APPENDIX I 


—-COOP »sA DORSEYs S/1S/25933500092. 
—-FTNsoLsE-« 


CGA OV VOLVO NAO OA OERADN ANNAN ANA ADADAANANADN ANNAN AANANAAN OM 


PROGRAM HRZGRD1 

THIS IS A PROGRAM TO COMPUTE THE BEARING AND DISTANCE 
ERRORS IN A SONAR FIX WHICH RESULT FROM HORIZONTAL 
TEMPERATURE GRADIENTS IN THE WATER. 

THIS PROGRAM IS WRITTEN IN FTN-63 FOR A CDC 1604 COMPUTER. 

THE Y AXIS IS AN YSOVELOGITY LINE aikibge SHES VeLOCRIiY Is VO. 

THE ORIGIN (X=0sY=0) IS A POINT MIDWAY BETWEEN THE TWO 
SENSORS AND IS. THESINT ERSEGRIGNSOF THE Y AXISIAND ALINE 
THROUGH THE TWO SENSORSe 

SENSOR A IS LOCATED AT X1lsY1 AND SENSOR B IS AT XX1sYYle 

A Lie BEARING SeW Iii THE PREEIiXeBRGm@AREeREWATIVE TO AN AXIS 
THROUGH THE TWO SENSORSe 

AllweBEARIINGS DESTGNAIMEDRGY “GRiaks LETTERSMARE RELATIVE TO 
THE Y AXISe 

ALL RANGES AND DISTANCES ARE IN METERS 

ANGLES IS A SUBROUTINE TO CALCULATE lic mROM O TO 360 
DEGREES USING TANGENTS. 

BRGA IS THE BEARING FROM A TO THE APPARENT TARGET e 

BRGB IS THE BEARING FROM B TO THE APPARENT TARGET. 

BRGG IS THE BEARING FROM THE@O@RIGIN T@ -HE ACTUAL TARGET. 

BRGP IS THE BEARING FROM THE ORIGIN TO THE APPARENT TARGET 

BRGT IS THE BEARING FROM THE APPARENT TO THE ACTUAL TARGET. 

ERROR IS THE DISTANCE BETWEEN APPARENT AND ACTUAL POSITION. 

RACT IS THE ACTUAL RANGE OF THE TARGET FROM THE ORIGIN. 

RAPP IS THE APPARENT RANGE OF THE TARGET FROM THE ORIGIN. 

RAPPINC IS THE RANGE INCREMENT. 

RAPPMAX IS THE MAXIMUM APPARENT RANGE CONSIDERED. 

REDRNG IS APPARENT RANGE DIVIDED BY SENSOR SEPARATION.’ 

RELERR IS THE RELATIVE DISTANCE ERROR (ERROR/RAPP) e 

RNGERR IS THE DIFFERENCE BETWEEN RACT AND RAPP. 

RPD CONVERTS ANGLES IN DEGREES TO ANGLES IN RADIANS. 
(RADIANS=DEGREES TIMES 201745329252 RADIANS PER DEGREE) 

5S IS SENSOR SEPARATION. 

VO IS THE VELOCITY AT X=0s VGRAD 1S THE VELOCITY GRADIFNT. 

XACTsYACT IS ACTUAL TARGET POSITION. 

XAPPsYAPP IS APPARENT TARGET POSITION. 

KINC IS THE DISTANCE INCREMENT GN THE Deeb PRE CTION. 

YDIFF IS THE DESIRED PRECISION OEN CLOSGieGiee OF HE 
COMPUTED AND ACTUAL POINTS OF. INTERSECTION OF TWO RAYS. 
ALPHA IS THE BEARING OF THE APPARENT TARGET FROM SENSOR Ae 
BETA IS THE BEARING OF THE APPARENT TARGET FROM SENSOR Be 

DELTA IS THE DIFFERENCE BETWEEN GAMMA AND PHI. 

GAMMA IS THE BEARING FROM THE ORIGIN TO THE ACTUAL TARGET. 

OMEGA IS THE ANGLE BETWEEN THE Y AXIS AND A LINE THROUGH 
THE TWO SENSORSe IT CAN VARY FROM O TO +90 DEGREES. 

PHI IS THE BEARING FROM THE ORIGIN TO THE APPARENT TARGET. 

PI IS 341415926536 

TAU IS THE BEARING FROM THE APPARENT TO THE ACTUAL TARGET. 
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READ IN DATAe 


100 RPD = .0)] 745s 252 $ P1=3e1415926536 

110 READ lils VOsVGRADsYDIFF 

111 FORMAT(3F10.0) | 

120 READ 121s SSsRAPPsRAPPMAX sRAPPINCsXINC 

121 FORMAT (5F10.0) 

PROGRAM STOPS WHEN SENSOR SEPARATION IS NEGATIVE. 

122 IF (SS) 900,900,125 

125 OMEGA=0.0 

132 OMEGAR=OMEGA#RPD 

140 Y1=SS/2*COSF (OMEGAR) $ YYl=-Y1l 

141 X1=SS/2*SINF(OMEGAR) $ XXl=-xl 

170 REDRNG=RAPP/SS 

183 PRINT 184sO0MEGA SS sRAPPsREDRNG, V0 sVGRADsYDIFF 

184 FORMAT(1H1s16(/)s5X»s8HOMEGA = F3e0s9X95HSS = F6e039Xy% 
L7THRAPP = F7.0s9X%s9HREDRNG = F4.39//15X%s5HVO = F6els 
212Xs8HVGRAD = F7.5912Xs8HYDIFF = F3els//) 

185 PRINT 186 

186 FORMAT (2Xs4HBRGP4X » 4HRACT 6X 9 2HV26X 9 4HBRGA4X » 4HBRGBSX 9 
1 5SHERROR4X » 4HBRGT 3X » SHDELTA3X » GHRNGERR 3X 9 6GHRELERR ) 

200 DO 600 I[=1535 ‘ 

201 IF(I-18) 205,520,205 

ROUTINE TO SET APPARENT TARGET POSITION 

205 BRGP. = 10.*! 

206 ICOUNT=0 y 

210 PHI=OMEGA+BRGP 

TEST TO ENSURE THAT PHI IS LESS THAN 360 DEG.~ 

211 IFCPHI~360.)21352133s212 

212 PHI=PHI-360. 

ROUTINE TO COMPUTE RELATIVE BEARING OF TARGET FROM SENSORS 

213 PHIR=PHI*RPD 


214 XAPP=RAPP*SINF(PHIR) S YAPP=RAPP*COSF (PHIR) 
220 K=1 & CALL ANGLES(X1sY19XAPP»YAPP sALPHA 9K) 
221 ALPHAR=ALPHA#RPD SBRGA=ALPHA-OMEGA 


222 \PRBRGA) 2 2eece occ) 

223 BRGA=BRGA+360. © 

225 K=2 $ CALL ANGLES (XX1sYY1sXAPPsYAPP »BETA 5K) 
226 BETAR=BETA*RPD S$BRGB=BETA-OMEGA 
227 IPC SRGS). 22052 505-0 

228 BRGB=BRGB+360.- 


ROUTINE TO TRACE RAYS AND COMPUTE ACTUAL TARGET POSITION 


230 XA=X1 $$ YA=Y1 . $ XB=XX1 &S YB=YY1 
231 IF(OMEGA) BB252993 232 
232 VA=V0+VGRAD#X] $ VB=V0+VGRAD*XX1 
ALPHAR1=ALPHAR S$ BETARI=BETAR 
TEST TO SEE IF X1 IS LESS THAN XAPP IN ABSOLUTE VALUE 
233 IF(ABSF(XAPP)-ABSF(X1)) 24092803280 
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ROUTINE TO COMPUTE XACTsYACT WHEN X1 IS GREATER THAN XAPP 
IN ABSOLUTE VALUE. 
240 RATIOA=COSF(ALPHAR1)/VA $ RATIOB=COSF(BETARI1)/VB 
X22XX1 S$ XXINC#XINC 
241 X2=X2+XXINC $ V2=VO+VGRAD*X2 


IFCABSF(V2*®RATIOA) eGTele) 2575242 
242 IFC(ABSF(V2*RATIOB) eGTele) 2575243 
243 ALPHAR2=2-*PI-ACOSF(V2*RATIOA) 
BETAR2Z=ACOSF (V2*RATIOB) 
DXA=XA-X2 $ DXB=X2-XB 
NOTE SIGNMOF DRA 15 AEWA +. 
DYA=DXA/TANF( (CALPHARI+ALPHAR2) /2.¢) 
DYB=DXB/TANF ((BETARI+3ETAR2)/2.) 
Y2A=YA-DYA Seer? B=WBT DVS 
TEST FOR INTMRSEGIION OF RAYS 
GHECK TO SEEWIF CRITIGAREL.ANGLE Ut Boilie. 
ICOUNT=ICOUNT+1 S$ IF(ICOUNT.LTe¢200) 2445530 
244 IF(180.-BRGA) 245352059255 
245 IFC Y2A-Y2B) 24652583241 
246 IF(C(Y2B-Y2A)eLTeYDIFF) 2585257 
255 IF(Y2B=Y2A) 25625 Cece 
256 IF((Y2A=Yaepmebi«YDIFE) 25652 ae | 
257 X2=X2-XXINC S$ XXINC=el¥*XXINC $ GO TO 241 
258 XACT=X2 $ YACT=Y2A ‘S$ GO TO 400 


ROUTINE TO GET XB EQUAL TO XA WHEN ABSOLUTE VALUE OF X1 IS 
LESS THAN ABSOLUTE VALUE OF XAPP 

280 IF(PHI-180.)281592813291 y, 
CALCULATION OF STARTING POINTS OF RAYS WHEN PHI IS LTe 180 
281 DX=XA-XB 

BETAR1=ACOSF(VA*COSF( BETAR)/VB) 

DY=DX/TANF( (BETAR+BETAR1)/2.) 
285 YB=YB+DY $ XB=XA $$ V1=VO+VGRAD*®XA 

ALPHA1=ALPHA $ BETAL=BETARL/RPD 
AT THIS POINT BOTH RAYS HAVE THE SAME X COORDINATE. 
286 GO TO 300 


CALCULATION OF STARTING POINTS OF RAYS WHEN PHI IS GTe 180 
291 DX=XA-XB 
NOTE THAT SIGN OF DX IS ALWAYS +e 
ALPHAR 1=2-*PI-ACOSF(VB*¥COSF(ALPHAR)/VA) 
DY=DX/ TANF ( (ALPHAR+ALPHARI1)/2.) 
295 YA=YA-DY $ XA=XB $ V1=V0+VGRAD*xXB 
ALPHA1=ALPHARI1/RPD $ BETAL=BETA 
AT THIS POINT BOTH RAYS HAVE THE SAME X COORDINATE. 
296 GO TO 300 


299 V1=VO $ ALPHAR1=ALPHAR $ BETAR1=BETAR 
VA=VO+VGRAD¥X1  $ VB=VO0+VGRAD¥XX1 
ALPHA1=ALPHA $ BETA1=BETA 


AT THIS POINT THE RAYS ARE DISTANCE STEPPED IN THE + OR —X 
DIRECTION (DEPENDING UPON PHI) UNTIL THEY INTERSFCT, 


5) 


wd 
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XVIT 
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300 DX=0.0 $ Y2A=YA $$ Y2B=YB 

301 RATIOA=COSF (ALPHARI) /V1 S$ RATIOB=COSF(BETAR1)/V1 
XXINC=X INC 

305 IF CPHI «Li. 1S@mmg06+536 


CALCULATION OF XACT sYACT WHEN PHI IS LESS THAN 180 DEGe 
306 DX=DX+XXINC 
307 V2=V1+VGRAD#*DX 
IF(ABSF(V2*RATIOA) e-GTele) 327,308 
308 IF(CABSF(V2*RATIOB) Glele) 3275309 
309 ALPHAR2=ACOSF(V2*RATIOA) $ RETAR2Z=ACOSF(V2*RATIOB) 
DYA=DX/TANF ( (ALPHARI+ALPHAR2)/2.) 
DYB=DX/TANF ((BETARI1+BETAR2)/2,.) 
313 Y2A=YA+DYA $§ Y2B=YB+DYB 
IF (BREA S25 5 ores STS 
TEST FOR INTERSECTION OF RAYS 
314 IF(180.-BRGA) 32535203315 
315 IF(Y2A-Y2B) 31633283306 
316 IF ((Y26ZAP LE OVO PRE). 328 5 oe 7 
325 IFCY2B-Y2A) 326593283306 
326 ITF((Y2A=Y2Z0.LT eYDVEF 1328 9327 
327 DX=DX=-XXINC SB XXINC=e1l*XXINC $ GO TO 306 
328 XACT=XA+DX $ YACT=Y2A ~ $ GO TO 400 


CALCULATION OF XACTsYACT WHEN PHI IS GREATER THAN 180 DEG. 

336 DX=DX+XXINC 

NOTE THAT SIGN OF DX IS “AbWArs ey 

337 V2=V1-VGRAD#DX 
IFCABSF(V2*RATIOA) eGTele) 3575338 

338 IFCABSFI(V2Z*ERATIOS) -GTiele) 35 fuse 

339 ALPHAR2=2.*PI-ACOSF(V2*RATIOA) 
BETAR2Z=2e¢*PI-ACOSF(V2*RATIOB) 
DYA=DOX/TANF ( (ALPHAR1+ALPHAR2)/2¢) 
DYB=DX/TANF((BETARI+BETAR2)/2,.) 

NOTE THAT SINCE SIGN OF OX WAS KEPT + THROUGHOUT >» SIGhigE 

DYA AND DYSB WILL BE THE OPPOSITE OFOWHAT THEY SHOULD ME 

343 Y2A=YA-DYA $ Y2B=YB-DYB 4 

Tet OR INTERSECTION OF RMAs 

344 IF(1806e¢-BRGA) 34555203355 

345 IF(Y2A-Y2B) 34653585336 

3 SSE PFU IY 2ZE-Y 2A RY OIE) S5e0e 57 

355 IF (Y2B=2A) 356 54583336 

3565 RX ( YAReeS. LT o YD er ee Sis 3 5./ = . 

357 DX=DX-XXINC B XXINC=el*XXINC S GO TO 32@e 


358 XACT=XB-DX $ YACT=Y2A $ GO TO 400 
ROUTINES FOR CALCULATION OF ERRORS 
400 XERR=XACT-XAPP $ YERR=YACT+YAPP 


405 ERROR=SORTP(KERR*XERR + YERR*YERR) 
410 RACT = SQRTF(XACT#XACT + YACT#YACT ) 
415 RNGERR = RAPP = RACT 


54 


APPENDIX I 


C SUBTRACT RNGERR FROM RAPP TO GET RACTe 
416 RELERR=ERROR/RAPP 

C Z IS A DUMMY VARIABLE TO ENABLE USE OF SUBROUTINE ANGLES. 
420 Z=0.0 

XVITI 421 K=3 $ CALL ANGLES(ZZ XACT sYACT sGAMMA 5K) 

422 DELTA = PHI — GAMMA 

C SUBTRACT ERRORS FROM PHI TO GET GAMMA. 
425 K=4 $ CALL ANGLES(ZsZsXERRsYERR»TAUsK) 
427 BRGT = TAU — OMEGA 
428 IF(BRGT)4293500,500 
429 BRGT=BRGT+360. 

€ 

€ 

S 


CALCULATIONS COMPLETE» PRINT RESULTS. 


500 PRINT 510s BRGP»sRACT »V2sBRGA  sBRGBsERROR »BRGT sDELTAs 
IRNGERRsRELERR 

510 FORMAT (2X9F3e¢092X9F7e094X 9 F641 93X92(F5e193X) oF 760. 
13XsF5el93BX9F4els2XoF8e093XoF6,4) 

515 GO TO 660 

520 PHI=PH]I+1056. & PRINT 5215 PHI 

521 FORMAT(F5e0924Xs28HIMPOSSIBLE FIX(BRGP=180 DEG)) 
GO TO 600 

530 PRINT 531s PHI 

531 FORMAT (F5e0924Xs22HCRITICAL’ ANGLE LIMITED) 


f 600 CONTINUE 
AC 
€ 
C 
E 
C 


XIX 


INCREMENT RAPP AND GET NEW RESULTS. 


610 RAPP=RAPP+RAPPINC / 
Hest TO.see IF RAPP HAS. EXCEEDED RAPPMAX.N IF YESs READ 
NEW DATAs IF NOTs GET RESULTS FOR NEW RAPP. 
620 IF(RAPP-RAPPMAX) 17031703120 
900 CONTINUE 
» & END 
© 
C SUBROUTINE TO CALCULATE ANGLES FROM 0 TO 360 DEGREES 
C USING TANGENTS. | 
C 


SUBROUTINE ANGLES (AsBsCsDsANGLE 5K) 
700 RPD = 201745329252 
POA ee GD=8) 71057055710 
705 IF(C-A) 70697383707 


706 ANGLE=270. $ GO TO 740 
707 ANGLE=90. $ GO TO 7.40 
710 ANGLER= ATANF((C-A)/(D-B) ) $ ANGLE=ANGLER/RPD 


felboeer (NCA OB los? 16s tong 
HY OmIF (E-A)/ 203 7505 722 

Ti? IF(CHA) 72297365740 

aC iG OB 07305 1385 136 


"720 ANGLE=ANGLE+360. $ GO TO 740 
722 ANGLE=ANGLE+180. $ GO TO 740 
730 ANGLE=180. $ GO TO 740 


736 ANGLE=06 $ GO TO 740 
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738 GO TO( 750s 76097405770) sK 


740 RETURN 
TO PRIN? Ou 


S°GO TO 740 


751 FORMAT(35Xs18HTARGET AT SENSOR A) 


760 PRINT 761 


$ GO TO 740 


761 FORMAT(35Xs18HTARGET AT SENSOR 8B) 


T70 PRIN yr 
771 FORMAT 


END 


-EXECUTE. 
1500 
100000 
200000 
400000 
~99%9 « 


©0001 
100000 
100000 
100000 
~999 6 


$ GO TO 740 
(35Xs19HTARGET AT XAPP,YAPP) 


al 
S000 00 1qn000 
1000000 1n0000 
2000000 200000 
a aed -~999 6 


19000 
19000 
10000 
=~99% « 


APPENDIX [1 


uN Read in Data:VO, VGRAD, YDIFF 
Read in Data 
SS, RAPP, RAPPMAX, RAPPINC, XINC 


No 
: 
Yes 


Eounure Sensor Position 
(2) 3 Print Input 2 Print Output -Headings 


Set ae Target Posit.— 
Compute Target Bearings 
Set Starting 


VIL —— 81 XIII} Coordinates 
and Velocity 










Yes = NO 
VIII tt 7 


ABSF(XAPP) -ABSF(X1) <0 





No Le Yes 
PHI > 180 


Calculate Y Coordinate 
of Ray B When X = XAI 

Step X Toward Xl “Until 
Rays Intersect at XACT, 
ACT 













IX 





Calculate Y Coordinate 


Acree VY COOEGimace XII 
of Ray B When X = Xl of Ray A When X = XXI 
XIv} set Starting ve lbeity 
and Coordinates 
No - Tf Yes 
PHI > 180 


Step X in -X 
XVII| Direction until Rays 
Intersect at XACT, YACT 







) (a) 











Step X in +X 
XVI Direction until Rays 
Intersect at XACT, YACT 


s| XVIIT Compute Errors 
Print Results 


If XACT, YACT have been 
Bound, ctOgeal ter aL 









Yes = 
Increment RAPP No 
. L£ RAPP > .RAPPMA ~ 


Dy . : 
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TABLE 4 


Effect of Sensor Separation (SS = 1 km.) 
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Effect of Sensor Separation (SS = 5 km.) 
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Effect of Sensor Separation (SS = 10 km. ) 
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TABLE 9 


Effect of Target Range (RAPP = 1 kn. ) 
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TABLE 10 


Effect of Target Range (RAPP = 3 kn.) 
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Effect of Target Range (RAPP = 100 km.) 
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TABLE 13 


Effect of Target Range (RAPP = 2,000 km.) 
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TABLE 14 


Effect of Sensor Axis Orientation (OMEGA = 0°) 
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(OMEGA = 10°) 
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